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Fermions, as a major class of quantum particles, provide platforms for quantum information 
processing beyond the possibilities of spins or bosons which have been studied more extensively. 
One particularly interesting model to study, in view of recent progress in manipulating ultracold 
fermion gases, is the fermionic version of measurement-based quantum computation (MBQC), which 
implements full quantum computation with only single site measurements on a proper fermionic 
many-body resource state. However, it is not known which fermionic states can be used as the 
resource states for MBQC and how to find them. In this paper, we generalize the framework 
of spin MBQC to fermions. In particular, we provide a general formalism to construct many- 
body entangled fermion resource states for MBQC based on the fermionic projected entangled 
pair state representation. We give a specific fermionic state which enables universal MBQC and 
demonstrate that the non-locality inherent in fermion systems can be properly taken care of with 
suitable measurement schemes. Such a framework opens up possibilities of finding MBQC resource 
states which can be more readily realized in the lab. 



I. INTRODUCTION 

Quantum computation can be realized with differ- 
ent quantum degrees of freedom, for example photons 
and spins. Fermions, as another major class of quan- 
tum particles, have been relatively less explored for 
their application in quantum computation and can lead 
to new possibilities. Although it is expected pQ that 
fermions have polynomially equivalent quantum compu- 
tation power as spins/bosons, it is possible that sub- 
exponential speedups can be achieved with fermions over 
spins/bosons in certain computational tasks. For ex- 
ample, the quantum simulation [2j [3] of fermionic many- 
body systems can be much more easily implemented with 
fermionic degrees of freedom due to the intrinsic sign is- 
sue in the simulation of fermions with spins/bosons. 

The possibility of using fermions for quantum com- 
putation has been studied in a few contexts. It has 
been shown that the circuit model quantum computa- 
tion can be implemented with fermions which efficiently 
simulates quantum circuits with spins [1. On the other 
hand, topological quantum computation can be real- 
ized using certain two dimensional fermion states with 
strong correlations^ j. In particular, it is known that 
the fractional quantum Hall state with filling fraction 
v = 5/2 can support universal topological quantum 
computation^. Moreover, quantum teleportation[6], an 
important quantum protocol for both quantum commu- 
nication and quantum computation, have also been gen- 
eralized to fermion systems [7]. 

The measurement-based quantum computation 
model [11 has been extensively studied in spin systems 
where quantum computation is implemented with 
only single spin measurements on a proper many- 
body entangled resource state. Many spin resource 
states are known [ITTH2Q] but a large scale experimental 
realization has not been achieved. With the recent 
exciting experimental progress in manipulating ultracold 
fermion gases [8- 10 , it is then interesting to ask whether 



similar computational schemes could be implemented in 
fermion systems, with only single site measurements on a 
fermionic resource state which ideally can be realized in a 
controlled way with ultracold fermionic atoms. With the 
large variety of quantum states that exist in simple free 
fermion systems, like Fermi liquids, quantum Hall states, 
and topological insulators, a fermionic version of MBQC 
may provide new platforms for quantum information 
processing with reduced experimental complexity while 
at the same time enjoying the same advantage as in the 
spin MBQC model that no coherent quantum operations 
are needed to carry out the whole computation. 

However, no theory exists for fermionic MBQC which 
studies what fermionic resource states are useful and 
what single site measurement patterns are necessary to 
achieve universal quantum computation. Naively, one 
might expect that a direct Jordan Wigner mapping of 
spin resource states to fermions would give a useful 
fermionic resource state for MBQC, but this is not true as 
the mapping is nonlocal and local spin measurements on 
the resource state can no longer be implemented with lo- 
cal fermion measurements after the mapping. Moreover, 
one of the key properties wanted for a MBQC resource 
state is lost during this mapping. It is highly desirable 
to have the MBQC resource states be the ground states 
of local Hamiltonians and many spin resource states are 
designed to have this property[15-20 . Unfortunately 
this property is not preserved by the nonlocal mapping 
to fermions and the resulting fermion states can no longer 
be generated by engineering the appropriate local Hamil- 
tonian terms in the system and then lowering the tem- 
perature. Therefore different approaches are needed to 
construct a useful fermionic MBQC model. 

In this paper, we show that MBQC is possible in lo- 
cal fermion systems by presenting an explicit construc- 
tion of a fermionic resource state together with the single 
site measurement patterns necessary to realize univer- 
sal quantum computation. Our construction is based on 
the fermionic Projected Entangled Pair States (fPEPS) 
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representation [2TJ [22], which is known to describe ground 
states in local fermion systems [23]. The construction gen- 
eralizes the idea of designing spin MBQC resource states 
based on the spin PEPS representation [T3] [T4] to fermion 
systems. By encoding the quantum information to be 
processed into the even parity sector of local fermion 
modes, we demonstrate how universal quantum compu- 
tation can be achieved on a fermionic state with only 
single site measurements. One complication arising from 
this encoding is the extra fermionic measurement pos- 
sibilities in the odd parity sector which introduces non- 
local by-products to the computation. We demonstrate 
further that such by-products in the computation can 
be properly taken care of by keeping a 'fermionic' frame 
of the by-products together with the Pauli frame as in 
the spin MBQC models [TTJ [12]. Starting from this ex- 
plicit construction, we expect that the fPEPS formalism 
could yield fermionic resource states with simpler encod- 
ing scheme and as the ground states of more easily real- 
izable local Hamiltonians. Therefore, we also discuss in 
general how to design fermion resource states for fMBQC 
from fPEPS representation. 

The paper is organized as follows: In section [II] we 
review how measurement-based quantum computation is 
done in spin systems. In particular, we focus on the in- 
terpretation of measurement operation on the resource 
state as teleportation in the virtual space in the PEPS 
representation, which is crucial for our generalization to 



fermion systems. In section III, we start from the basic 
building block of MBQC-teleportation, and show how it 
can be realized in fermion systems. Putting the telepor- 
tation steps together, we obtain a simple fermionic re- 
source state in section [IV] A- B and demonstrate in detail 
how each step in MBQC can be realized on such a state. 
This example is the starting point of a more general con- 



discuss following Ref. [13] how such teleportation steps 
occur in the virtual space of the PEPS representation 
when measurement operations are performed on individ- 
ual spins in a many-body entangled resource state. This 
line of thought allows us to generalize the MBQC scheme 
from spins to fermions. 



A. Teleportation 

Teleportation is a way of achieving quantum compu- 
tation by doing multi-spin measurements and it can be 
thought of as the basic building block for MBQC which 
realizes the universal quantum computation with single 
spin measurements. Originally, teleportation was discov- 
ered as a way of transmitting information [6 - that when a 
measurement is done in the Bell basis, information stored 
at one place can flow to another. Later on, it was found 
that not only the information but also some extra gates 
can be teleported as well if measurements are done in a 
modified Bell basis [26]. 

First, to see how teleportation transmits informa- 
tion, consider the following setup that involves an in- 
put qubit to be teleported and an entangled pair 
\E) = 1 00) + 1 11) (suppressing normalization) shared be- 
tween the input end and the output end. The information 
in is transmitted to the output end when one mea- 
sures and half of the entangled pair jointly. Further- 
more, by choosing different measurement bases, different 
gates can be performed on while it is teleported to 
the output end. 

To see this, consider an input qubit 



l^i) =m |0)+mi|l), 



(1) 



struction based on fPEPS which we present in sectionjTV] with |m | 2 + |mi| 2 = 1, and an entangled pair | £^23)) 
" ' 1 shown in Fig. [l] Qubit 1 and 2 belong to the input end. 

Qubit 3 belongs to the output end. 



C. Finally, we conclude and discuss future directions in 
section W\ 



II. REVIEW: MEASUREMENT-BASED 
QUANTUM COMPUTATION 

In this section, we review the measurement-based 
quantum computation (MBQC) scheme, which claims 
that universal quantum computation can be achieved 
with only single-site measurement operations on a many- 
body entangled spin state. We focus on the interpreta- 
tion by Gross and Eisert[13] of the computational power 
of resource states in terms of their Projected Entangled 
Pair State (PEPS) represent at ion [21 ] [25 ] , In such a rep- 
resentation, measurements on the physical lattice sites 
in the resource state correspond to teleportation steps in 
the virtual space of entangled pairs, which can then be 
composed and designed to simulate full quantum circuits. 

Following this logic, we first review how teleportation 
with entangled pairs can perform a universal set of uni- 
tary gates while transmitting the information. We then 




FIG. 1. Teleportation of one-qubit unitary gates. The qubits 
are denoted by dots and the entangled pair is depicted by 
a dashed line. The big circle represents the input end and 
qubits inside it are measured together. After measuring qubit 
1,2 together information in qubit 1 flows to qubit 3. 

The total wave function of the system is 

fe} = (m |0) + mi |1» (|00) + |11)). (2) 

If we measure qubits 1 and 2 in the Bell basis |</>i 2 ) = 
a a (8) 7(|00) + |11)), a = 0,1,2,3, where a a are Pauli 
matrices, the wave function of the unmeasured qubit 3 
results in (^l^m) = o- a (mo |0) + mi |1)). Thus, one 
can see that the original information is teleported from 
qubit 1 to qubit 3 with possible extra Pauli operations. 
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In general, teleportation not only transmits informa- 
tion, but also implements gates at the same time, as 
the Pauli gates seen at the output in the example given 
above. Measuring in a basis of the generic form 

|0? 2 ) = (^V a )0/(|OO) + |ll)), (3) 

where U is any one-qubit unitary gate, yields at the out- 
put 

(^2^123) = a a U(m |0) + mj |1» 3 = a a U . (4) 

For example, we can choose U as the Hadamard opera- 
tion H = |0) (+| + 11) (-1, where |±) = ^(|0) db |1», and 

the phase gate Z{6) = e _i i |0) (0| + e l i |1) (1| to imple- 
ment the corresponding one-qubit gates during telepor- 
tation. 

In addition to one-qubit gates, the requirement of uni- 
versal quantum computation also involves certain two- 
qubit gates, for example, the controlled- Z gate together 
with the Hadamard gates on the two qubits [27 which we 
denote as U p h- Specifically, 

U ph = |00) (++| + |01) (+-| + |10) (-+I-I11) (—1 . (5) 

The schematic for teleporting U p h is depicted in Fig. |2j 
It can be interpreted as a generalized version of teleporta- 
tion which involves a two-qubit input and three entangled 
pairs [28], It can be checked that measuring qubits 1,2,3 
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FIG. 2. Teleportation of the two-qubit controlled operation 
Uph- As shown in the figure, qubit 1 is the control qubit and 
qubit 5 is the target qubit. After proper 3-qubit measure- 
ments are implemented on qubits 1,2,3 and 4,5,6 respec- 
tively, Uph gets teleported from qubits 1 and 5 to qubits 7 
and 8. 

in 

(a x y®(<T X y ®i(\o + +)±\i--)), (6) 

and measuring qubits 4,5,6 in 

(a x ) k ® (a x ) 1 ® I(\00+) ±\11-)), (7) 

with k,l = 0,1, teleports U p h to qubits 7,8 up to 
Pauli operations on 7,8 separately. Hereby, one can see 
that teleportation is indeed a way of realizing universal 
quantum computation, with multi-qubit measurements. 



B. Projected Entangled Pair States and spin 
Measurement Based Quantum Computation 

A projected entangled pair state (PEPS) [25] is a way 
of expressing a many-body entangled state as a projec- 
tion from a product of maximally entangled pairs. PEPS 
was first invented to study many-body systems in con- 
densed matter theory; however, it turned out to be use- 
ful for understanding the power of resource states in 
measurement-based quantum computation. More explic- 
itly, if we imagine the maximally entangled pairs in PEPS 
as in a virtual space where teleportation can be achieved, 
and interpret the physical Hilbert space of spins as a pro- 
jection from many virtual spins, then measurements on 
the physical spins correspond to teleportation steps in 
the virtual space and it may be possible to implement 
a universal set of unitary operations on virtual qubits 
by merely performing single-qubit measurements in the 
physical space [T3]. 

For simplicity, let us first consider PEPS on a one- 
dimensional chain. As depicted in Fig. |3j a spatially 
one-dimensional (ID) virtual space is a chain consisting 
of maximally entangled pairs shared between nearest- 
neighbor sites. With D-dimensional virtual spins, the 
maximally entangled pairs are in state ^f^ 1 At 
the left and right end of the chain, there are boundary 
states \L) and \R). On every site, there are two virtual 
spins, each being half of an entangled pair connecting 
neighboring sites. Shortly, we discuss how virtual and 
physical spins are related via a projection on each site. 
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FIG. 3. Illustration of a ID PEPS. The virtual space con- 
sists of left and right boundaries |L), \R), and virtual spins 
entangled in Xw_o The big circle represents an on-site 

projection where virtual spins inside are projected together 
to the physical space. 

The wave function for the virtual chain follows 

wv)=\l) n ■ (t>> M+1 V> (8) 

k=l \i=0 J 

where k labels different sites in the chain. 

A PEPS with (i-dimensional physical spins is obtained 
by a local projection P on virtual spins located on each 
site in the virtual space which maps the virtual space to 
physical space. In Fig. [3J the projections are presented 
as circles. Note that the projection here is only a map 
between two Hilbert spaces and is not the usual sense of 
projection that needs to obey P 2 = P. More specifically, 

N N-l /D-l \ 

\4>peps) = n ^ w n e mm+i \ r ) ( 9 ) 

k=l k=l \ t=0 ) 
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with local projection operators on each site defined as 

d-l D-l 



ft = EE \?) A Ur(er\. 

i=0 t,r=0 



(10) 



Note that each i) is a state in the d-dimensional phys- 
ical space and A\ r ^ are coefficients of P that depend on 
|i) in the physical space that one wants to project onto. 
Also note that the physical dimension refers to the inter- 
nal degrees of freedom of spins and is different from the 
spatial dimension of the lattice. 

Using Eq. [9]&[l0j one can show that the wave function 
of a ID PEPS with N sites can be expressed as 

d-i 

\^peps) = < L I A n \ R ) 1^2 In) • 

i k =0 

(ii) 

The PEPS construction provides a perspective to see 
the relation between MBQC and teleportation [13]. Here 
we give an explicit example that illustrates this idea. 
Consider measuring site 1 in \<j)) = a 5) + b |l) (d=2). 
After the measurement, the wave function of the unmea- 
sured physical spins becomes 



(0\^PEPS) = 



¥ A l )Ai 2 . 



•Aft\R)\i 2 h iiv). 



The form of the state remains unchanged while the left 
boundary gets teleported to site 2 and is changed into 
(L\ (a* A + b* A 1 ), which can be viewed as a gate act- 
ing on the old boundary. From Eq. 10 , we learned that 



the physical and virtual space are related by the projec- 
tion. Given a certain P, we can relate measurements in 
the physical space to teleportation in the virtual space as 
well as changes in lattice boundary to unitary operations 
teleported in the virtual space. Quantum computation 
could therefore be achieved in the virtual space by choos- 
ing appropriate measurement bases in the physical space 
that correspond to bases for teleporting a universal set 
of gates in the virtual space. 

In the above, we have seen that unitary operations on 
one spin could be realized with MBQC on a ID PEPS. 
With suitably chosen projection P, arbitrary single spin 
operations could be implemented. Yet, for the universal- 
ity of quantum computation, entangling operations be- 
tween two spins must also be feasible; thus a more gen- 
eral 2D lattice is required for MBQC. For this purpose, 
a 2D PEPS as shown in Fig. [4] can be constructed simi- 
larly. The only differences are that the virtual space now 
contains both vertical (between sites and (i + 1, j)) 
and horizontal (between (i, j) and (i, j + 1)) entangled 
pairs ^SzQ 1 N)' an d the on site projection P^ becomes 



d-l D-l 




pa = E 



E 

u,l,r,d= 



\hj) A \j,uird {ulrd\ 



(12) 



FIG. 4. Representation of a 2D PEPS with input modes on 
the left boundaries. 



which is almost the same as the ID projection in Eq. 
[T0| except that the number of virtual indices is doubled. 
Note that i denotes the state of a physical spin as defined 
in Eq. [TO} 

One of the best known universal resource states for 
MBQC is the 2D cluster state, with on site projection 

Puidr = |6) (00 + +| + |l) (11 1. It can be checked 

that measurement on this projected space can be used to 
input information, implement the universal set of gates 
discussed in section |II A[ correct possible computation 
by-products due to measurement randomness and finally 
read out the computational result [28], Therefore, the 2D 
cluster state can serve as a resource state for MBQC. Var- 
ious other resource states can be constructed using this 
framework [T3l fT4| [29] . but a general rule is still missing 
for determining whether a projected entangled pair state 
can be used as a resource state or not. 



III. FERMIONIC TELEPORTATION 

Many-body fermion systems exist naturally in avail- 
able experimental settings and it would be nice to gener- 
alize the mechanism of MBQC to fermions. Our goal in 
the next sections is to show that on a many-body fermion 
state single site measurements can be used to simulate 
the result of any quantum circuit. We achieve this with 
a similar procedure as that used in spin MBQC: we first 
construct fermionic teleportation steps for implement- 
ing a universal set of gates and then use the fermionic 
PEPS (fPEPS) representation to map teleportation in 
the virtual space to single-site measurements in a physi- 
cal state. However, fermions are very different from spins 
in two specific ways: 1. fermion operators anti-commute 
with each other and fermion wave functions are anti- 
symmetric; 2. the total parity of a fermionic system is 
always preserved. Therefore, care must be taken in map- 
ping from spins to fermions and the generalization is far 
from direct. In this section, we start with a fermionic 
version of teleportation, discuss the necessity of a new 
encoding scheme for fermionic systems, and finally give 
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a way to achieve universal quantum computation with 
fermionic teleportation. These serve as the basic building 
blocks for fermionic MBQC discussed in the next section. 



A. Fermionic teleportation as a generalization of 
spin teleportation 



carried throughout the whole scheme 2. information is 
encoded in a fixed parity sector and all operations have 
fixed parities. In the following, we apply these strategies 
to the general cases. First we need to specify the en- 
coding scheme of general quantum circuits into fermion 
modes. 



In quantum computation with spins, the analogy to 
the bits 0/1 in classical computation is the two- level spin 
up/down states |0) / |1), or qubits. In the fermionic case, 
the information is encoded in the wave function of local 
fermionic modes. It seems straight forward to define the 
two-level states by the occupation number of the modes. 
Namely, the analogy of |0) is a state with no fermion in a 
mode, or vacuum and that of 1 1) is a) |fi), where a) is 
the creation operator for a fermion mode. The maximally 
entangled spin state 1 00) + |11) can be replaced accord- 
ingly by two entangled modes defined as (1 + aja^) |^)- 

However, this naive mapping fails as one attempts to 
do fermionic teleportation with the configuration shown 
in Fig. [I] where each dot now represents a fermion mode. 
As fermion parity of a system is always preserved, the 
input mode 1 cannot be in a superposition state of (mo + 
mioJ) In order to deal with this problem, we take 
a route similar to that in Ref. [T| and [7| by adding extra 
modes and encoding information in a fixed-parity sector. 
As depicted in Fig. |5j we add an extra mode and a second 
pair of entangled modes, such that the input defined in 
Eq. [l] becomes | ^13) = (mo \Q) + mia\a\ |^))i3, which 
has a definite even parity. To check the feasibility of 
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FIG. 5. Teleportation of two-mode fermionic unitary gates. 
Fermionic modes are depicted as white dots, and entangled 
mode pairs are represented by a dashed line. Input modes 
(1, 3) are in state 1^13) = (wio |^) +m\a\a\ where mode 
1 carries the information, and mode 3 preserves the parity of 
the mode 1. 

this strategy, we show in the following that the input in 
modes 1 and 3 can be teleported to modes 5 and 6: The 
total wave function of the system is 

1 

= J2^a(44) a (al4 + l){a\a\ + 1) . (13) 

a=0 

A measurement of modes 1-4 in (</>| = (Q\ (1 -\-a4a3a2a1) 
results in state (<p\ip) = (mo + m\a\a\) \Q) on mode 5 
and 6, which is exactly what is desired. 

This is the simplest case of a teleportation circuit and 
illustrates the general strategy we take to deal with the 
special property of fermions: 1. a proper ordering of all 
fermion modes needs to be given at the beginning and 



B. n — ► 2n encoding scheme 

Due to the parity constraint discussed above, extra 
modes are needed when encoding spin states into 
fermion modes to preserve the total fermion parity 
of a system. Various encoding schemes have been 
proposed [1 which satisfy this constraint. For discussion 
in this paper, we choose to encode 1 qubit into 2 
fermionic modes, or more generally, n qubits into 2n 
fermionic modes. As illustrated in Fig. [5] and Fig. 
|6j in our scheme, a 'parity' mode is assigned to every 
'info' mode containing the real information to ensure 
that the total parity of an info mode and the auxiliary 
parity mode is always fixed. Thus, a spin system with 

n qubits in state \i/j n ) = J2{ ai } m {a,} \ a i a 2 a n ), 

where {^} = {ai,a2, ,& n } is encoded into 

a fermionic state with 2n modes \^ n ) f = 

E M m {ai M a l P ) ai (44 P ) a2 -4aWu P r-\^}, 

where i p is the parity mode of the info mode i. Note 
that here we have chosen the order of the modes such 
that fermionic operators a\ always appear in front of 
otj for i < j, and the fermion parity of the state shall 
always be even. 

As a spin state with n qubits is encoded into a fixed 
parity fermionic state with 2n modes, spin gates must 
also be redesigned accordingly so that an n-qubit spin 
operator is encoded into a 2n-mode parity preserving 
fermion operator and the universal set of spin gates are 
mapped to a set of fermionic gates which possess the 
same universality. 

A generic one-qubit unitary spin operator 

U = Uoo |0> <0| + U w |1> <0| + U01 |0> <1| + U n |1> <1| 

(14) 

is encoded into a 2-mode fermionic gate (where mode 1 
is the 'info' mode, and mode 2 is the 'parity' mode of 
mode 1): 

Uf = L^00^1^I^2^2 + ^01 tt l^2 + ^10 a l 

With this encoding, we have the 2-mode fermionic 
phase gate 

Zf{&) = aia\a2al + e %e a\aiala2, (15) 
and the fermionic Hadamard gate 

Hf = aia{a2<^2 + a \ a 2 + a i a 2 — q4<^ic4<^2- (16) 
which can be composed to simulate arbitrary unitary 
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gates on a single qubit. 



Fermionic teleportation for a universal set of 
gates 



Similarly, the 2-qubit U p h is mapped to U p hj on 4 
consecutive fermionic modes, modes 1, lp, 2, 2p, where 
1,2 are control and target modes, and lp,2p are parity 
modes of 1,2 respectively. For simplicity, here we de- 
note fermionic states a 1 and (1 + ajaj) as 
|0}j,|l)p and (|00) + |11))/ where we always order the 
modes as 1, lp, 2, 2p, etc. 



U ph .j = |0000) / «00| 
+ |0011) / «00| 
+ |1100) / ((00| 
-11111)^(001 



<n|)/«oo| 
<n|)/«oo| 

(H|)/((00| 
(H|)/((00| 



(111)/ 
(111)/ 
(111)/ 
(111)/- 



(17) 



Therefore, with Eqs. [l5j [16] and [lT| a universal set 
of fermionic gates for fermionic quantum computation is 
constructed and can be used to simulate the universal 
set of spin gates for the original spin quantum computa- 
tion. Note that the gates discussed here only act on the 
even fermion parity sector and are unitary only within 
this sector. However, as information is encoded fully in 
this sector, these gates are sufficient for quantum com- 
putation and are implemented in the MBQC scheme de- 
scribed below. Unlike in the fermionic circuit model of 
quantum computation [1 where fully unitary fermionic 
gates are necessary, in fermionic MBQC simulating such 
quasi-unitary operations is sufficient and can be readily 
realized. The odd fermion parity sector contributes to the 
fermionic MBQC scheme as computational by-products 
when the measurement result falls into this sector. As 
we show below, such computational by-products can be 
properly dealt without destroying the universality of the 
computation scheme. 



Now that we have defined the encoding of states and 
the mapping between gates, in the following, we show 
that the universal set of fermionic gates can be imple- 
mented by measuring entangled fermionic states in cer- 
tain bases, thus achieving universal quantum computa- 
tion with fermionic teleportation. Note that in the dis- 
cussion of this section, we always assume that each pair 
of 'info' mode and 'parity' mode always have even par- 
ity. The occurrence of odd parity pairs is considered as 
computational by-products later. 

The schematic for teleporting an arbitrary two-mode 
parity preserving fermionic gate [//, the equivalent of a 
1-qubit spin gate, is shown in Fig. [5] By comparing Figs. 
[T] and [5j we can see that the number of inputs as well as 
that of entangled pairs are both doubled in the fermionic 
case. The wave function of the state that corresponds 
to Eq. [2] in the spin case is given in Eq. [l3j It can be 
checked that the measurement on modes 1-4 in basis 
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Uoia\a\ 



• Uioa\al + Uiia\alala\) \Q) 

(18) 

teleports Uf to mode 5 and 6. 

As for teleporting the 4-mode controlled operation 
Uph-j, the setup depicted in Fig. [6] is utilized. The wave 




HO o- 

12 



IB 
-O 

■O 

16 



FIG. 6. Teleportation of the 4-mode controlled operation. In- 
put modes 1 and 9 are the fermionic analogues of the control 
and target, and modes 3 and 11 are the parity modes corre- 
sponding to 1,9 respectively. After measurements on modes 
1^6 and 7 ^ 12 are implemented, the gate is teleported to 
modes 13 ~ 16. 

function of this state is 



m a6 (aI4) a (4«I 1 ) 6 (l + 4a\ 3 )(l + a\a\ A ){\ + a\a\){\ + a\a\){\ + a\ a\ 5 )(l + al 2 a f 16 ) . (19) 



a, 6=0 



One can check that U p h-j can be teleported by measuring modes 1-6 of the top site in 
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)* = (!- «2« 



4 



«5 a 6 



t^t 



ai a. 



1"3 



(20) 



and mode 7-12 of the bottom site in 



= (1" 



-a\ a\ 2 - 



- QlyQigQlQ 



' (21) 



Hereby, we have successfully found a measurement 
bases corresponding to a universal set of gates and 
have shown that universal quantum computation can be 
achieved by teleportation with fermions. 

However, our consideration so far is over simplified as 
we have assumed that the computation always occurs in 
the even fermion parity sector and the measurements al- 
ways result in the basis we want. In fact, measurement 
errors always occur as we can not choose which particular 
basis among a complete set to measure in. Measurement 
errors in teleportation steps lead to unwanted by-product 
operations being teleported. For fermion states, it is also 
possible to change the parity sector of the states, which 
seems to pose a serious problem for our scheme. We ad- 
dress these issues in the following sections and show that 
they can be properly taken care of and will not impede 
our ability to do MBQC. We refer to the extra operations 
teleported as 'by-products' instead of 'errors' to empha- 
size that the former is due to the intrinsic randomness of 
quantum mechanics and cannot be avoided while the lat- 
ter is due to noise and perturbation and can in principle 
be reduced. 



A. A simple example of fermionic resource state 

Here we demonstrate that fMBQC can be achieved us- 
ing a special fermionic resource state on the lattice shown 
in Fig. [7[ Like what we have seen in fermionic teleporta- 
tion, the number of input and entangled pairs are doubled 
in the fermionic case compared to the spin case; thus, the 
spin lattice for MBQC (shown in Fig. [4|, which has 3 or 4 
qubits on every site corresponds to the fermion lattice in 
Fig. [7] which has 6 or 8 modes per site and two entangled 
pairs connecting neighboring sites. 



(1,2) 

(id ° 

{lf±} o o- ■ o o- 

?? ?? 



(2,1) 



(4,1) 



<>-.<M>-. 



o 66 o- ■ o 66 o- 





FIG. 7. Representation of the 2D fPEPS simple example re- 
source state. The lattice consists of boundary modes on the 
left and entangled mode pairs connecting every site. 



IV. FERMIONIC PROJECTED ENTANGLED 
PAIR STATES FOR MBQC 



Even though we have demonstrated the viability of 
fermionic teleportation for individual gates in the pre- 
vious section, our ultimate goal is to show that a 
circuit consisting of multiple operations can be simu- 
lated with local measurements, or in other words, to 
achieve fermionic measurement-based quantum compu- 
tation (fMBQC). Thus, it is necessary to have a fermionic 
lattice state similar to the spin lattice state in Fig. [4] 
which allows multiple steps of measurements as the in- 
formation flows from one place to another. 

In this section, we first assemble the teleportation steps 
and give a simple yet universal example of resource state 
for fMBQC. We examine in detail the possible by-product 
operations that occur in the measurement process and 
show how they can be taken care of with proper mea- 
surement schemes. We then discuss the more general 
fermionic Projected Entangled Pair States (fPEPS) for- 
malism which, like PEPS for spin, allows more possibili- 
ties for finding novel fermionic resource states. 



The lattice consists of input mode pairs aa's on the 
left boundary and entangled pairs (/?J -at \Q} and 

{f3'\ -a' J 7+1 + l) \ for horizontal bonds, and (tI^J+i 7 + 



1)|0) and (YJ ; S 



,t r/t 

,3° 



1) |0) for vertical bonds. The 
labeling of modes on a site is shown in Fig. [8] When 
writing the bonds, we define the ordering of sites on the 
lattice as left to right (i, j + 1) and top to down 
(i + 1, j). This state can be thought of as a fermionic 
PEPS with a trivial projection on each site. In the fol- 
lowing we think of all the modes on each site as one big 
degree of freedom and discuss how MBQC can be imple- 
mented with single site measurements on this state. 

To see the feasibility of fMBQC in this example, we 
give in detail the procedure to implement each necessary 
step in fMBQC on this state. 

• Assume WLOG that the input modes on the left 
boundary are all initialized in |fi), and measure- 
ments are performed on the sites in the first column 
from top to bottom, and then column by column 
from left to right so that the information flows to 
the right. 
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FIG. 8. Labeling of modes on a site. Modes on the left are 
labeled a,a\ modes on the right are f3/3\ modes at the top 
are 55' , and modes at the bottom are 77'. 



• Just like in MBQC for spins, the lattice is initially 
entirely entangled. As one wants to achieve cer- 
tain operations, for example, two-mode gates or 
the four- mode controlled operation U p h-j, which in- 
volve only one or two entangled rows, one needs 
to isolate the rows and decouple them from other 
rows. To isolate a row, we remove its entanglement 
with neighboring upper and lower rows by measur- 
ing modes on the sites of the upper and lower rows 
in the occupation number basis 



't' 
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for all n = 0, 1. Apply this measurement wherever 
necessary to prepare the lattice for the implemen- 
tation of a particular circuit. 

• After partially decoupling the lattice when neces- 
sary in the way introduced above and using the 
results from Eqs. [THJ [20j [21] we see that we can 
implement a universal set of fermion gates by sin- 
gle site measurements on this state. But this is 
not enough to claim universality for MBQC as we 
have not considered the effect of measuring in basis 
other than the desired one. We discuss how to deal 
with the computational by-products introduced by 
the randomness in fermionic measurements in the 
next section. 

• We finally read out the output on the right bound- 
ary by measuring the sites in the occupation num- 
ber bases. Therefore for fermions, we just mea- 
sure the rightmost column in (a^a^) Ua \Q) or 



a 1 a 
suits. 



^ Ua |^), with n a = 0, 1 to yield the re- 



B. Dealing with measurement randomness in the 
simple model 

In this section, we address the effect of measurement 
randomness in our scheme. As discussed in the previous 
section, fMBQC could in principle be achieved with mea- 
surements in certain bases; however, measurement results 
in orthogonal bases that span the rest of the Hilbert space 
may lead to by-products to the simulated operation. This 



can be viewed as the fermionic analog of the Pauli by- 
products that emerge when a Bell measurement is per- 
formed. In general, we cannot choose which basis state 
results from the measurement and whenever a measure- 
ment is done, a by-product occurs. Therefore, dealing 
with by-products becomes a necessity to make sure that 
there is a finite probability of simulating the wanted oper- 
ation in order to achieve efficient quantum computation. 

So far in our discussion, we have used two important 
assumptions: 1. we required that the pair of input modes 
on a site always have even parity as we designed one mode 
as the parity mode of another; 2. we only mentioned the 
measurement basis that gives rise to the desired answer 
without discussing other orthogonal bases that would po- 
tentially produce by-products. In general, the fermion 
parity constraint only requires that the total parity of 
modes be fixed. Therefore, the input mode pairs could 
also have odd parity, eg. m^a^ + mia^ Similarly, 
there are no other constraints on the measurement bases 
as long as the total parity is fixed. As a result, it may 
seem that the choice of bases is arbitrary, leading to all 
kinds of by-products in the simulated operation. Yet, 
for the consistency of our scheme, we choose a complete 
set of measurement bases in which the mode pairs (ac/, 
/?/?', 77', 55') each have a fixed parity, which could be 
either odd or even. Therefore, depending on the par- 
ity of the measurement bases, we could characterize the 
by-products into two categories: 

• Parity-preserving by-products: by-products that 
come from measurements which preserve the parity 
of the input (i.e. with an even parity measurement 
basis), such that the parity of the output is the 
same as the input. These parity preserving local 
fermionic operations can be mapped to local spin 
operations and could be corrected locally as spins. 
To see this more explicitly, we assume on a ID chain 
with input of even parity, measurement in the state 
(1 + aa'Pfi') \Q) simulates the desired operation. 
Other orthogonal bases that preserve the parity are 
(l-aa'P?) (aa'+pp') and (aa'-pp') \Q). 
It is obvious that the by-products on the output for 
these bases are the fermionic equivalent of Pauli Z, 
X, and Y respectively. In our simple example, such 
by-products can all be incorporated into the next 
operation to be teleported and hence get corrected. 

• Parity- violating by-products: by-products that 
emerge in measurements which change the parity 
of the input (i.e. with an odd measurement basis), 
such that the parity of the output is the opposite of 
the input. Using the example above, the orthogo- 
nal bases in the odd sector are (aa'P ± p') \Q) and 
{aa! ft' d= P) To keep the information flow, a 
corresponding encoding of spin states into the odd 
parity sector needs to be defined, for example by 
requiring that the 'parity' mode always has the op- 
posite parity to that of the 'info' mode. Moreover, 
this type of by-products are not the typical spin by- 
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products. Instead they implement odd fermionic 
operations on the input, which maps back to non- 
local spin operations. Nevertheless, since in our 
scheme we have required a fixed parity on each 
mode pair, the nonlocal part of an odd operation 
only contributes an overall (±1) to the total state, 
and therefore we only need to worry about the local 
part which is correctable by local measurements. 
Note that this is a special property of this exam- 
ple. Generally by-products from odd measurement 
bases are non-local, and we discuss the general case 
in the next section. 



C. General fPEPS construction for fMBQC 

1. Review of f PEPS formalism 

In the previous section, we demonstrated that fMBQC 
is feasible in principle on a 2D lattice. However, in this 
model, the on-site measurements involve many degrees of 
freedom and the resource state may not be readily realiz- 
able. Our ultimate goal is to find a resource state which 
contains few modes per site and is the unique gapped 
ground state of a simple Hamiltonian, for example a free 
fermion Hamiltonian. The computational power of such 
a state is connected to its physical properties through 
its fPEPS representation. fPEPS, like PEPS for spin, 
represents many-body fermion states as projections from 
entangled virtual fermion pairs and provides new possi- 
bilities for finding fermionic resource states for quantum 
computation. 

First, we review the fPEPS formalism [2T1 [22]. A 2D 
fPEPS is obtained from a lattice of fermionic entangled 
pairs (for example as shown in Fig. [9] or Fig. [7]) by pro- 
jecting the fermion modes on each site to a smaller phys- 
ical Hilbert space. For the simple example given above, 
the projection is trivial on each site. In a general fPEPS 
state, the boundary modes and the entangled modes be- 
tween sites (a, /?, 7, 5) are only virtual and we denote the 
physical modes as c to distinguish them from the virtual 
ones. The virtual entangled mode pairs are again ordered 
from left to right (i, j + 1) and top to bottom 
(i+l, j). The virtual boundaries and mode pairs between 
sites are denoted as B iy i and H\- = {ol\^0\^ +1 + l)fc \Q) 

for horizontal bonds and Vjj = (7} jSj+i j + l)fc |^) for 
vertical bonds respectively, where the integer k labels the 
number of bonds per direction per site. Fig. [7] and Fig. 
[9] represent models with k = 2 and k = 1 respectively. 

The wave function for the virtual space can be ex- 
pressed as [21] 

|V>„ = II B >i"u V U \ n )v ■ (22) 

An on-site projection P^ that maps the virtual space 
to the physical space with physical modes q on every site 



(U) (1,3) (1,4) 
(1,1) o o- ■ -o o- ■ -o o- ■ -o 

? ? ? 

SV\ SV\ 

(2,1) O O- ■ -O O- - O O- ■ -O 
? ? ? 

ri\ /^\ /"~\ 

(3,1) O O- • -O O- • O O- ■ -O 
? ? ? 

(4,1) O O- • -O O- • -O O- ■ -O 

FIG. 9. Representation of a 2D fPEPS with k = 1. The 
lattice has only one bond per direction. 

is defined as [21]: 

Pij = £ ^[w]n(4"c^r< Qfc 7r).. 

{n}=0 l,k %3 

(23) 

where {n} is the set of occupation numbers for every 
mode, and A^[{n}] depends on the intrinsic properties 
of the physical state one wants to project to. 

Pij is constrained to have a fixed parity for the result- 
ing state to be physical, or 

^2(ri£ + np k + n lk + n ak + n 5k ) mod 2 = c, (24) 

l,k 

where c is constant for each site. 

To yield a physical state, one applies the projection 
operator P^s to the virtual state \i(j) v together with the 
physical vacuum state and then takes the vacuum 
expectation value on the virtual space as all the virtual 
modes must be annihilated and only physical modes are 
left, 

W P = v(^\l[p ij \^ v \n) p 

= v <fi| II P n II B '^ Il U V U \ Q )v P)p (25) 

This form of many-body fermion state is the starting 
point for a more general construction of fermionic re- 
source states for MBQC. 



2. Information flow in fPEPS 

In the following, we show how the information stored 
in the left boundaries B^i gets transmitted to the right 
by local measurements on an fPEPS. For simplicity, from 
now on, we assume k = 1, as shown in Fig. [9j The mea- 
surements are performed site by site from top to bottom 
and then from left to right starting from site (1,1) in 
column 1. 

To illustrate the flow of information, we first look at 



10 



the measurement on site (1,1). Suppose site (1,1) is 
measured in = On |^n) p . To see what the state 

becomes after the measurement, we first rearrange \ip) p 
in Eq. [25] and commute the terms containing modes on 
site (1,1) together. We get: 

a (i,j)^(l,l) 

Q a ltl = J2 v(^u\PiiBl 1 H 11 Vu\^ii) v \n 11 ) p (26) 



eral formulation, it is possible to find physically more fea- 
sible fPEPS resource state for MBQC. Extra care needs 
to be taken when dealing with measurement randomness 
on a general fPEPS and we discuss briefly possible diffi- 
culties in the next section. 



3. Dealing with measurement randomness for general 
fPEPS 



where a denotes different terms in Bi^ if input modes on 
site (1,1) is entangled with other modes of B^i, which 
could possibly happen in a generic state as long as the to- 
tal parity is fixed. Note that no extra signs are produced 
in this procedure as we are free to move the entangled 
pairs and the projections because they have fixed pari- 
ties. 



Using Eq. [25] and Eq. [26] it is obvious to see that after 
the measurement, the state becomes 

l^ll) <01l|^) p = x 

v (^1 ll(i,j)^(l,l) P ij R l,l B i,l H ij V ij \ty v \®)p > 

with 

R1i= P <fin I «<nn|Oiii'iiBi iHiiVii|nii)jnii> 1> . 

(27) 

R% 1 is an operator on the a virtual mode of site (1,2) 
and the 5 virtual mode of site (2, 1). We can hence in- 
terpret the effect of measuring site (1,1) as information 
flow in the virtual space from site (1,1) to site (1,2) and 
(2,1) as shown in Fig. 10 The encoded state changes 



from Bi^i to R\^\ and correspondingly certain operation 
is implemented. This is similar to the picture we had 
with spin MBQC where measurements on physical sites 
correspond to operations implemented on the informa- 
tion flow in the virtual space. 



(1,2) (13) (1,4) 
(1,1 o o- -o o- -o o- -o 



(U) (1,3) 



(3,1 



o /o 
(2 l o o— ■ o 




o— ■ o o— ■ o o- ■ o 

{3/1 |o^.(^.(oV. 
(4,1 



FIG. 10. Illustration of boundary changes after measurement 
on site (1,1). The boundary originally on site (1,1) moves to 
site (1, 2) and site (2, 1). The new boundaries are Bi^(i / 1) 
and Rn. 

This formalism provides a general framework to study 
MBQC based on many-body fermionic state. The simple 
example we studied before falls into this framework with 
k = 2 and trivial projection = /. Based on the gen- 



In dealing with measurement randomness for the sim- 
ple model, we classified the by-products into two cat- 
egories depending on whether the output contains the 
same parity as the input. Yet, we showed that the by- 
products are local and locally correctable as all the mode 
pairs (aa', /?/?' ', etc) have a fixed parity. However, in 
a general fPEPS state with possibly an odd number of 
bonds per direction {k = 2n+l), n G N) and a more gen- 
eral encoding scheme, non-local by-products could occur 
and special attention is needed when designing MBQC 
schemes based on such states. As we show in the follow- 
ing, the non-locality of such by-products can be properly 
taken care of with careful design and is not a fundamen- 
tal difficulty in using fPEPS states as MBQC resource 
states. 

We use the k = 1 model to illustrate the basic idea. 
Assume that we are measuring the sites in a column- wise 
order, i.e. we first measure the first column from first 
row to last row and then second column from first row to 
last row, etc. Let us look more closely at the measure- 
ments on column 1, starting from site (1,1), and moving 
downward. Suppose that the boundary modes are always 
ordered from up to down. So after measurement on site 
(1, 1), they are ordered from site (1, 2) to site (2, 1) to site 
(3, 1) etc. The parity constraint with a general encoding 
scheme is that the whole boundary chain has a fixed total 
parity, but each boundary mode may not. In particular, 
the boundary mode on site (1,2) might not have a fixed 
parity. This leads to extra sign effect when site (2,1) is 
measured. In particular, if site (2, 1) is measured in an 
odd basis which corresponds to an odd operation on the 
boundary, it applies a non-trivial sign factor (— 1) Uq i,2 to 
the boundary mode on site (1,2). Similarly, measuring 
site in an odd basis causes a non-trivial sign fac- 

tor on sites (i',2) for i' < i. Therefore, the by-product 
induced is indeed non-local. 

In general, after finishing measurements on the j'th 
column, the overall sign accumulated on site 

+ 1) in column j + 1 is determined by the number of 
odd measurement bases below site (i, j) in column j and 
the occupation number operator n Q; ( ij+1 ) of a mode on 
site (i, j + 1). Define 



N, 



(28) 



i' ,i'>i 



where is the parity of the measurement basis on site 
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Then, we obtain 

S (iJ+1) = (29) 

Even though the by-products are non-local, they can 
be dealt with in a local way. Note that as long as one 
keeps track of all the measurement results, the total by- 
products that happen to the boundary modes can be de- 
termined after one finishes the measurements of one col- 
umn. Moreover, the by-products factorize into a product 
form, of individual operators on each boundary mode 
separately, for example as given in Eq. |29j Such by- 
products can be incorporated into the operation to be 
implemented when measuring the column j + 1 and can 
be corrected locally just like correcting Z by-product in 
spin systems. 

To sum up, in a general fMBQC scheme based on 
fPEPS, non-local by-products do occur. But as they fac- 
torizes into a product form, they can be corrected locally. 

V. CONCLUSION 

In this paper we generalized the measurement-based 
quantum computation scheme from spin systems to 



fermion systems. We gave a simple example of many- 
body fermion states and demonstrated how it could be 
used as a universal resource state for MBQC. More gener- 
ally, we provided a framework for constructing fermionic 
resource states for MBQC based on the fPEPS repre- 
sentation of fermion states and discussed ways to deal 
with the non-local by-products that might come up in 
the general scheme. 



This framework provides a general starting point for 
the construction of new MBQC schemes. The ultimate 
goal is to find resource states that are easy to realize 
experimentally, for example in a free fermion system 
where particles move around but do not interact with 
each other. Unlike spin systems, which factorize into 
total product states and lose all computational power 
without interaction, the hopping of fermions in the lat- 
tice and their non-trivial mutual statistics can generate 
entanglement among different sites in space, which can 
subsequently provide the basis for the power of quantum 
computation. 
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